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摘 要 : 利用 广义 p 值 和 广义 置信 区 域 的 概念 对 含有 三 个 随机 效应 的 Panel 数 据 模 型 中 回归 系数 的 假设 检 
验 问题 建立 了 精确 检验 ， 构 造 了 回归 系数 的 几 个 广义 置信 区 域 ， 讨 论 了 本 文 所 构造 的 检验 和 置信 
区 域 在 尺度 变换 下 的 不 变性 ， 对 这 几 种 检验 的 功效 和 置信 区 域 的 覆盖 率 给 出 了 数值 模拟 结果 。 
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Panel 数 据 模 型 本 质 上 是 一 类 有 具有 套 误 差 结构 (nested error structure) 的 线性 回归 模型 ， 
广泛 应 用 于 计量 经 济 学 等 领域 ， 该 类 模型 的 研究 不 仅 受 到 经 济 学 家 也 受到 了 统计 学 家 的 广泛 
关注 ， 近 年 来 有 不 少 这 方面 的 研究 成 果 发 表 。Panel 数 据 模型 回归 系数 的 广义 最 小 二 乘 F 检 
验 (文献 中 记 作 Fgzrs(p)) 中 含有 未 知 参数 ，Rao 等 由 对 参数 进行 估计 并 给 出 一 种 可 行 性 FF 检 
验 ， 有 些 文献 中 则 忽略 群 内 相关 性 采用 标准 的 最 小 二 乘 下 检验 F Wl MFHT BIER 
出 一 种 新 的 检验 Fa4(p)。 这 些 方法 都 采用 了 下 检验 ， 但 都 不 是 精确 检验 ，Wang 和 Liskilal 对 它 
们 功效 分 别 进行 了 研究 。Wang 和 Ma 还 对 此 模型 提出 了 两 个 精确 的 下 检验 Fp 和 Fyw， 并 利 
用 Fs 和 Fw 构造 出 一 个 检验 。 

本 文 将 基于 广义 p 值 检验 理论 方法 给 出 回归 系数 的 几 个 新 的 精确 检验 和 置信 区 域 ， 并 通过 
对 功效 和 置信 区 域 履 盖 率 的 模拟 对 它们 进行 比较 。 广 义 p 值 检验 最 初 是 为 了 解决 单 边 假设 检验 
中 由 于 元 余 参 数 存 在 而 无 法 获得 精确 检验 的 问题 而 提出 的 。Tsui 和 Weerahandil®) 和 Weeraha- 
ndil6] 分 别提 出 了 广义 p 值 (generalized p value) 和 广义 置信 区 间 (generalized confidence in- 
terval) 的 概念 。 实 际 应 用 表明 在 精确 下 检验 不 存在 或 元 余 参 数 存 在 时 ， 基 于 广义 p 值 和 
广义 置信 区 间 来 获得 精确 检验 和 置信 区 间 的 方法 是 富有 成 效 的 ， 文 献 中 已 有 不 少 研究 成 
R14), Weerahandi”! 分 别 对 单 向 分 类 模型 方差 分 量 的 单 边 假设 和 多 向 分 类 模型 中 方差 分 量 
的 单 边 假 设 问题 给 出 了 基于 广义 p 值 的 精确 检验 。Zhou 和 Mathewl8] 把 这 种 方法 应 用 到 混合 效 
应 模型 中 ， 分 别 对 单个 方差 分 量 的 显著 性 和 和 两 个 独立 平衡 模型 方差 分 量 的 比较 建立 了 精确 
检验 ， 并 把 部 分 结果 推广 到 非 平 衡 情况 。Weerahandi fll Berger"), Chi 和 Weerahandilol 利 
用 广义 p 值 对 协 方差 具有 组 内 相关 结构 的 简单 生长 曲线 模型 中 回归 系数 建立 了 精确 检 
验 ，Lin Al Lee!) 将 他 们 的 结果 推广 到 具有 等 相关 结构 的 简单 生长 曲线 模型 。Ye Al Wangi? 对 
平衡 数据 下 一 般 随 机 效应 模型 的 单个 方差 分 量 构造 出 精确 检验 和 置信 区 间 ， 并 将 结论 推广 
至 两 个 独立 平衡 模型 方差 分 量 的 比较 。Tsui 和 Weerahandig 还 讨论 了 广义 p 值 检验 的 不 变 
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性 。Weerahandil6,13] 构造 了 不 同 指数 族 均 值 差异 性 的 广义 置信 区 间 ， 并 移 造 了 混合 模型 方差 
分 量 差异 的 广义 置信 区 间 。Krishnamoorthy 和 Mathew 利用 广义 枢 轴 量 对 对 数 正 态 分 布 的 
均值 建立 了 置信 区 间 。 有 一 些 文献 45,16] 是 对 广义 p 值 方法 进行 探讨 。 


2 广义 p 值 和 广义 置信 区 间 


设 2 为 分 布依 赖 于 参数 (7, 6) 的 随机 变量 ， 其 中 为 检验 参数 。56 为 元 余 参 数 ， 且 5 可 以 为 
参数 向 量 。 假 设 要 检验 的 问题 是 


Ho:n<m +— Hi:7n> m, (1) 


这 里 mo 为 预先 给 定 的 值 。 记 zz 为 Z 的 观测 值 。 
定义 1 设 全 (2,z,m,6) 为 随机 变量 Z，Z 的 观测 值 z， 及 参数 (7, 6) HER. FT 满足 
1) 1T(Z,z, 70,6) 的 分 布 与 见 余 参数 无 关 ; 
2) Ty(Z, 2,7, 6) 的 观测 值 T(z,z,m,6) 与 未 知 参数 无 关 ; 
3) MRE 25, T1(Z,z,n, 6) 的 分 布 关 于 7 随机 单调 增 或 随机 单调 减 ， 
WET (Z, z,m;,6) 为 一 个 广义 检验 变量 。 
1 (Z, z,n, 6) 的 分 布 关于 7 随机 单调 增 时 ， 对 于 检验 问题 (1) 定义 广义 p 值 为 


p = Pr(T%(Z, z,n, 6) > Ti(z,z,n, 5) | n = n0). 
471(Z, 2,7, 6) 的 分 布 关 于 了 随机 单调 减 时 ， 对 于 检验 问题 (1) 定义 广义 p 值 为 
p = Pr(T (Z, z,,6) < Ti(z, z,n, 6) | n = no). 


对 于 给 定 水 平 a， 如 果 p > a 则 接受 原 假设 ， 反 之 拒绝 。 由 条 件 1),2) 可 知 ， 广 义 p 值 与 任何 未 
知 参数 无 关 ， 所 以 是 可 以 计算 的 。 

EX2 设 宛 (2,z,7 5) 为 随机 变量 Z，2Z 的 观测 值 >， 及 参数 (7, 6) 的 函数 。 若 罗 满足 

1) ”了 T2(2,z,n,6) 的 分 布 与 未 知 参数 无 关 ; 

2) Te(Z,z,n, 8) 的 观测 值 T(z,z,n,6) 与 元 余 参 数 无 关 ， 
则 称 22(2, >z,m, 6) 为 一 个 广义 枢 轴 量 。 

对 于 感 兴趣 参数 的 广义 置信 区 间 可 以 利用 72(2,z,7, 6) 来 构造 。 例 如 ， 当 22(z,z 6) = 
7 时 ， 车 T(1 一 a) 表示 D2(2,z,n,6) 的 1 一 a 分 位 数 ， 则 ZT2(1 一 a) 就 成 为 了 的 广义 置信 上 限 。 
广义 置信 下 限 和 双边 置信 限 类 似 可 得 。 详 细 讨 论 参见 Tusi 和 Weerahandils]，Weerahandil6]。 


3 ”模型 分 析 
本 文 考虑 含有 三 个 随机 效应 的 Panel 数据 模型 
Yit = bo + wiBit++ cite Be + pit Acten, t=1,---,N; t=1,---,T. (2) 
这 里 ys 表示 第 i 个 个 体 在 时 刻 t 因 变量 的 观测 值 ，zitj 表示 第 ;个 个 体 上 第 5 个 自 变量 在 时 


At HUH: 6 = (B1,… ,Bk) 为 未 知 回归 系数 ， jii 为 第 i 个 个 体 效应 ， 假 定 这 NN 个 个 体 
是 从 一 个 大 的 总 体 中 随机 抽取 的 ，ji 是 随机 的 ; 和 表示 时 刻 t 的 时 间 效 应 ， 也 视 为 随机 
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的 ;es 为 随机 误差 。 假 定 所 有 pi, de 和 si 都 彼此 独立 ， 且 让 ~ N(0,07), 和 ~ N(O, 0), Eit ~ 
N(0,02)。(2) TEA 


其 中 


则 有 


其 中 


y = 1NnTho+ X8 +u, (3) 
y= (yn, try YT, Y21 YT, YN1 t’ YNT) ， 
X = (ai, ++, EIT, Par, ++ Bary ++ BNI» TNT), 
Tie = (vier, Diez, t, Lith), 
u= (In @17)p+ (An 8 Ir)à +£; p= (m, NY, 
A= (àn, Ar)’; e = (el1) +, err, END i , ENT) . 


Ds Cov(u) = o? (In Q Jr) + o3 (Jy Q Ir) + o2Inr, 


Jp =1717, JIn=1nl1ly, Inr =In @Ir. 


Jp =Jr/T, In =JIn/N, En =In- Jn, Er= lIr- Jr, 
Cov(u) = 02(Ey ® Er) + (To? + 02) (En 8 Jr) + (Nos + 02) (Un 8 Er) 
+ (No? +02) (Jn ® Er) 


E02Q1 + 02Q2 + 03Q3 + 07Q4, (4) 


oi = 02， o3 = To? + 02, o = No? +02, of = To? + Noi +07, 


Qi1= En ®Er, Q2=En®Jr, Q3=Jn Er, Qa = Jn 8 Jr- 


引 理 1 1) Qi, Q2, Q3, Qa We SH 

2) Qi, Q2, Q3, Q4 两 两 相互 正 交 ; 

3) Q1,Q2,Q3, Qa 的 秩 分 别 为 (N 一 1)(T 一 1), N-1, T-1,1. 

该 引 理 可 依据 王 松 桂 等 i" 中 相关 结论 加 以 证 明 。 分 别 用 Q, Qz, Q 左 乘 原 模型 ， 可 得 


yı = Qiy = Qı Xb + Qiu ê QIXB + u, 
Y2 = Qoy = Q2X P + Qou = Qi XB + ua, (5) 
yı = Q3y = Q3XB + Qau = Q3XB + u3, 


其 中 Eui = QiEu = 0, Cov(u;) = 02Qi， i = 1,2,30 假定 (X'Q;X)71, i= 1,2,3 均 存在 ， 这 在 
实际 数据 观测 中 是 容易 满足 的 。(5) 各 子 模型 中 可 得 到 6 的 最 小 二 乘 估 计 为 


Be = (X'QiX) 1 X'Qiy, i=1,2,3. 
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由 最 小 二 乘 统 一 理论 知 Bi 分 别 为 上 述 三 个 模型 中 6 的 最 佳 线 性 无 偏 估计 ， 且 
Cov(Bi) = 0?(X'Q:X) 1 £ 3i(o?), i=1,2,3. 
Fl GRAM a; = Qiy -—QiX Bi, BB) 0? 的 无 偏 估计 


üti y' (Qi — QiX(XQiX) IX'Qi)Y 


22 
6; = S? = 
ni Ni 


(6) 
其 中 ni = rk(Qi) —k, AY 


HAERES 中 相关 结论 可 得 。 

5/322 Âi, Bo, Bs, S2, S2, S2 相互 独立 。 

下 面 将 利用 本 节 获 得 的 优良 估计 来 构造 回归 系数 的 精确 检验 和 置信 区 域 。 
4 回归 系数 的 精确 检验 和 置信 区 间 


对 于 回归 系数 的 齐 次 假设 检验 问题 ， 本 节 将 提出 几 个 基于 广义 p 值 方法 的 新 的 精确 检验 ， 
并 且 这 些 检验 还 具有 良好 的 统计 性 质 。 考 虑 假设 问题 


Ho : HB =0 — H,: HB £0, (7) 
KH A Em x 大 矩阵， 其 秩 为 中。 沿用 上 一 节 的 估计 及 记号 ， 由 前 面 的 分 析 易 得 出 
ĝi ~ N(B, Di(o?)), Eio?) = 0? (X'QiX) +, HB. ~ N(HG,H(Ex(o7))H'), 1 =1,2,3. 
Ws? WA S?,i1=1,2,3 Ww, EX 


T; = (H,)' [H(£:(0?)) H] (HÊ:) — (Hbi)' [H (Si(o?s?.S7?))H"| (Hbi), 


= Ai — (Hbi) [H(2i(nis?V,-1))H"] (Hb), i = 1,2,3. (8) 


KH, b = (bi ,bik) HE AMM. BH p 

1) HB =O, A = (HB) [H (£:(0?))H' HÂ) ~ X2» Âi Vi = HF ~ ,相互 独 
立 ， 故 工 的 分 布 与 见 余 参数 无 关 ; 

2) 工 的 观测 值 t; = 0 与 未 知 参数 无 关 ; 

3) 又 工 为 HB; 的 正定 二 次 型 ， 故 其 在 备 择 假设 下 的 取 值 要 大 于 原 假设 下 取 值 。 
因此 ，(8) 对 假设 检验 问题 (7) 定义 了 三 个 广义 检验 变量 ， 据 此 得 到 广义 p 值 


pi = Pr(T, >0| HB =0) 


Pr(Ai > (Hbi) [H (Si(nis?V,-1))H’]’ (Hb:)) 


1 — By, {Fee [(H0;) [H(2i(nis7V,-')) A’) (Ab;)]}, (9) 


ll 


式 中 Fy: (-) 表示 X2 的 分 布 函数 ，Ev; 表 对 V 求 期 望 ，i = 1,2,3. 
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考虑 尺度 变换 
(8, ai) 二 (aß, aci), (Êi, Si) (aĝi, aS:), (a > 0). (10) 


可 以 看 出 ， 在 尺度 变换 (10) 下 检验 问题 (7) KAMET E b, S? 的 广义 检验 变量 工 均 保持 
不 变 ， 故 尺度 变换 (10) 下 假设 检验 问题 (7) 基于 广义 检验 变量 Ti = 1,2,3 的 广义 p 值 检验 是 
不 变 检 验 。 由 于 pi (i = 1,2,3) 定义 的 广义 p 值 为 度量 假设 问题 (7) 中 原 假设 正确 与 否 的 标准 ， 
即 如 果 广 义 p 值 小 于 预先 给 定 的 水 平 ag， 则 基于 广义 p 值 的 检验 方法 就 会 拒绝 原 假设 ， 反 之 接 
受 。 如 果 p; < wa， 这 说 明 从 现 有 数据 来 看 不 能 接受 原 假设 ， 此 时 自然 希望 能 够 构造 出 五 8 的 不 
变 置信 域 。 进 一 步 定 义 


T= (HB; — HB) [H (£:(s7?S2))H' (HÊ: — Hp), i=1,2,3. (11) 
由 于 
(HB; — HB) [H(2.(s;?S?))H’] (HÊ: — HB) 
= (HB; — HPY [H(£:(0?))H'] (HB; ~ HP). nis? V7 ', 
其 中 


x ` 2 
(HB; ~ HGY [H(£:(0?)) H] (HÂ; — HB) ~ Xh, Vi = TE o Aa, 
且 相互 独立 ， 可 知名 的 分 布 与 未 知 参数 无 关 ， 而 名 的 观测 什 
i = (Hb; — HPY [H(X'Q: X) H'] (Hb: — HP) 


与 见 余 参 数 无 关 ， 故 (11) 定义 了 三 个 广义 枢 轴 量 ， 并 且 在 尺度 变换 (10) 下 保持 不 变 ， 可 由 此 来 
构造 6 的 广义 不 变 置信 域 。 
对 于 非 齐 次 线性 假设 检验 问题 


Ho: H8 =d — Hi:HB#d 
这 里 及 定义 同 (8)， 且 HB = d 相 容 。 由 于 
E(HB;-d)=HB-d, Cov(HB;—d)= Cov(HB), i= 1,2,3. 
利用 (8) 的 有 关 结 果 ， 可 以 构造 广义 检验 变量 
T? = (HB; — d) [H (2:(0?))H'] (HB — d) — (Hb; — d) [H(2:(02s297?))H'] (Hb; — d). 
进而 得 到 相应 的 广义 p 值 计算 式 
pi = 1 - Ev{Fxs [(Hb; — d) [H(Zi(nis?V;-!))H’] (Hb; ~ qd)]}, i=1,2,3. 


其 中 各 符号 含义 同 (9)。 这 里 对 于 太 8 置信 域 可 同样 由 (11) 获得 。 
有 很 多 实际 问题 可 以 归结 为 单个 回 系 数 的 假设 检验 问题 ， 下 面 考虑 单个 同 归 系 数 的 单 边 假 
设 问题 的 精确 检验 。 记 By) 为 8 的 的 第 j 个 分 量 j e {1,2,.… ,k}， 考 虑 假设 问题 


Ho : By) < BO) — Hi : Bo) > Bey (12) 
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其 中 60) 为 一 预先 给 定 值 。 
对 于 假设 检验 问题 (12)， 定 义 


afo By) bia- Bo) 


Tio = 
o7 ))jj (2; (0?s?S;-7)) 55 


= Zi- Pau i= 1,2,3, (13) 

(Di(nas? Ve )); 

其 中 Bicy) ERABI ADE, bio) 为 Bigy WMA, (Di(o?)) 55 表示 Dilo?) 的 第 (j, j) TH 
Zi = (Big — Buy) (Bilo?));;) ~ N(0,1), Vi~, i=1,2,3 


且 相 互 独立 。 可 验证 Ti(j) 为 广义 检验 变量 。 用 B(-) 表示 标准 正 态 分 布 的 分 布 函数 ， 得 广义 p 值 
计算 式 


Pig) = Pr(Tio) > tio) | Bo) =L) 


Pr(Z > 0 ) 
(Zi(nis?Vi-*)) 55 


(14) 


| 


bis) — Be) }}. 
V (Silni? V); 
在 尺度 变换 (10) F, BRET RUE Â, S? 的 广义 检验 变量 Ti(y) 保持 不 变 ， 但 对 假设 检验 问 


题 (12) 本 身 的 参数 空间 发 生 了 变化 ， 所 以 假设 检验 问题 本 身 就 不 是 尺度 变换 下 的 不 变 检验 问 
题 。 考 虑 与 之 等 价 的 假设 检验 问题 


1- Ey, {| 


bo) Bip 0 
Ho: (5) = Fe < 9°) = ra «— Hy: bej) > 85): (15) 


对 假设 检验 问题 (15)， 定 义 


Pics 
TO i = 1,2,3. (16) 
(Dil(o?s? 15 


容易 验证 对 假设 检验 问题 (15)， 完 (j,i 二 1,2,3 为 广义 检验 变量 。 得 到 广义 p 值 计算 公式 
Buy = Pr(Tia > tw | O10) Nese) 
eo 
SERAI rar rc 外 
对 于 相应 的 尺度 变换 


(goi) — (gaoi)， (Bigg), Si) 一 (apiaSi)， (a> 0). 
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在 此 尺度 变换 下 ， 检 验 问 题 (15) 和 广义 检验 变量 (16) 均 保持 不 变 ， 因 而 基于 0;) 的 广义 p 值 
检验 为 尺度 变换 下 假设 检验 问题 (15) 的 不 变 检验 。 为 了 构造 关于 bo 的 不 变 置信 区 间 ， 考 虑 尺 
度 不 变 参 数 bij) = Bg s1 的 不 变 置信 区 间 。 定义 


Â, (q262972)).. 
(2 eg) (Pi(a; 87S; Djs 


Tx. 
+0) sı (2i(a2))77 


= Ziy (Di(nis?si2V 1D)77 ， 


RE ZEM, B Z ~ N(0,1)。 容 易 验 证 IT%,) 为 广义 枢 轴 量 ， 且 在 尺度 变换 下 保持 不 变 。 于 是 
可 以 利用 全; 的 分 位 数 来 构造 参数 0) 的 不 变 置信 区 间 。 


5 ”模拟 结果 和 模拟 方法 


本 节 首 先 给 出 对 假设 问题 Bo : HB =d — H : HBA d 基 于 T*, i = 1,2,3 的 广义 p 值 
检验 功效 的 模拟 结果 ， 记 8 = (H8 -d(H 一 d) 作 为 HB 与 4 之 间 偏 离 程度 的 一 个 度量 。 
表 二 给 出 了 基于 ,i = 12,3 构 造 瑟 6 的 广义 置信 域 对 其 真 值 覆 盖 率 的 模拟 结果 。 模 拟 中 
WH =(1,1,1), 8 = (12,1)， 设 计 阵 和 随机 产生 。 

由 表 1 的 结果 可 以 看 出 ， 本 文 构造 的 三 个 检验 的 功效 均 为 (H8 — d) (HB - d) 的 增 函 数 ， 其 
中 第 一 个 在 控制 犯 第 一 类 错误 概率 和 功效 增长 方面 都 明显 优 于 另外 两 个 。 表 2 可 看 出 本 文 所 构 
造 的 置信 区 域 基 本 上 包含 了 大 部 分 的 真实 参数 ， 但 从 稳定 性 及 和 置信 水 平 相 比 较 来 看 第 一 个 优 
于 其 它 两 个 。 因 此 ， 本 文 在 获得 的 三 个 精确 检验 和 置信 区 域 的 基础 上 ， 通 过 模拟 分 别 推荐 使 用 
第 一 个 检验 和 第 一 个 置信 区 域 。 


Al: 基于 TY 的 广义 p 值 检验 功效 (N =6, T=5, 0? =2, 02 =1, a=0.05) 


on tests b=0 b = 0.25 b = 0.64 =1 b = 1.44 b=4 
0.1 Tt 0.0590 0.1500 0.5900 0.8680 0.9840 1 
T 0.0300 0.0340 0.0500 0.0700 0.0880 0.1400 
T3 0.0400 0.0440 0.0480 0.0520 0.0500 0.1100 
0.2 T? 0.0590 0.1600 0.6020 0.8880 0.9800 1 
Ty 0.0220 0.0300 0.0560 0.0600 0.0700 0.1180 
T3 0.0460 0.0420 0.0440 0.0460 0.0600 0.1070 
0.5 Tř 0.0480 0.1580 0.6200 0.8820 0.9800 1 
TS 0.0320 0.0340 0.0500 0.0500 0.0460 0.1060 
T3 0.0480 0.0520 0.0540 0.0540 0.0660 0.1063 
1 Tt 0.0500 0.1550 0.6000 0.8800 0.9862 1 
T3 0.0213 0.0262 0.0275 0.0400 0.0450 0.1063 


T3 0.0488 0.0550 0.0563 0.0600 0.0612 0.2620 
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Ad: 互 B 的 置信 区 域 覆 盖 率 (N =T =5, o? =02 =1, 1-a=0.90) 


o2=1 o2=2 o2=3 o2=4 o2=5 o2 =6 
Ti 0.9455 0.9089 0.9616 0.9063 0.9505 0.8209 
Te 0.4036 0.8326 0.9674 0.1011 0.3424 0.7564 
Ts 0.6400 0.8811 0.8032 0.9116 0.2400 0.6621 


上 述 结果 的 具体 模拟 方法 

表 1: 

第 一 步 随机 产生 1000 个 Tx 的 观测 值 ， 统 计 其 中 大 于 0 的 比例 ， 此 比例 即 为 广义 p 值 ; 

第 二 步 重 复 第 一 步 1000 次 ， 对 给 定 水 平 a 统 计 小 于 a 的 p 值 比例 ， 此 比例 为 拒绝 原 假 设 的 


比例 ， 即 功效 。 


#2: 
IAT 的 分 布 产生 1000S OMT, MASEL, FE hE Ta 51-24 


位 数 之 间 ， 记 由 天 构造 的 广义 置信 域 覆盖 真 值 一 次 ; 


第 二 步 重复 第 一 步 1000 次 ， 统 计 和 覆盖 比例 。 
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Generalized p Value Test for Regression Coefficients 
in Panel Data Model 


CHENG Jing??, WANG Song-gui?, YUE Rong-xian’ 


(1- College of Mathematics and Sciences, Shanghai Normal University, Shanghai 200234; 
2- Department of Mathematics, Chaohu College, Chaohu 238000; 
3- College of Applied Sciences, Beijing University of Technology, Beijing 100022) 


Abstract: Some new exact tests and confidence regions of regression coefficients in the Panel data 
model with three random effects are established by using the concepts of generalized p value and 
generalized confidence interval. Invariance of these exact tests and confidence regions under the scale 
transformation are also discussed in the paper. The powers of these tests and coverage probabilities of 
these confidence regions are obtained through numerical simulations. 

Keywords: Panel data model; random effect; generalized p value; generalized confidence region; 
regression coefficient 


